First-principles study of the dipole layer formation at metal-organic interfaces 
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We study the dipole layer formed at metal-organic interfaces by means of first-principles calcula- 
tions. Interface dipoles are monitored by calculating the work function change of Au, Ag, Al, Mg 
and Ca surfaces upon adsorption of a monolayer of PTCDA (3,4,9, 10-perylene-tetra-carboxylic-di- 
anhydride), perylene or benzene molecules. Adsorption of PTCDA leads to pinning of the work 
function for a range of metal substrates. It gives interface dipoles that compensate for the differ- 
ence in the clean metal work functions, leading to a nearly constant work function. In contrast, 
adsorption of benzene always results in a decrease of the work function, which is relatively constant 
for all metal substrates. Both effects are found in perylene, where adsorption on low work function 
metals gives work function pinning, whereas adsorption on high work function metals gives work 
function lowering. The work function changes upon adsorption are analyzed and interpreted in 
terms of two competing effects. If the molecule and substrate interact weakly, the molecule pushes 
electrons into the surface, which lowers the work function. If the metal work function is sufficiently 
low with respect to the unoccupied states of the molecule, electrons are donated into these states, 
which increases the binding and the work function. 



PACS numbers: 73.30.+y,73.61.Ph,68.43.-h 
I. INTRODUCTION 

Applications of organic semiconduct ors in light- 
emitting diodeSjIi^ field-effect transistord^ and solar 
cell^^ have stimulated research into the fundamental 
electronic properties of organi c m aterials and their in- 
terfaces with metal electrodes .^^ISl The weak forces be- 
tween the molecules in an organic material lead to small 
band widths, which enhances the importa nce of elec tron- 
phonon and electron-electron interactions-I^ESlIIIill Nev- 
ertheless high charge carrier mobilities can be achieved 
in well-ordered molecular crystals.!^ As the quality of 
molecular crystals increases, transport of charge carriers 
across the interfaces between metal electrodes and the 
organic material starts to determine the performance of 
the devices.'i^ Metal organic interfaces (MOIs) often give 
rise to a non-Ohmic behavior, indicating the existence of 
significant Schottky barriers. 

Chemical bonding between molecules and metal sur- 
faces modifies the charge distribution at a MOI. It results 
in an interface dipole lay er, whic h strongly influences the 
Schottky barrier height .ISElIill This effect of chemical 
bonding is observed very clearly in self-assembled mono- 
layers of thiol ate molecules c hemically bonded to noble 
metal surfaces! ^^ * ^^ * ^^ * ^^ * ^° * ^^ * Common organic semicon- 
ductors however consist of closed shell molecules, which 
are usually thought to bind weakly to metal surfaces. 
It has therefore been assumed for a long time that the 
charge reordering at such MOIs is insignificant and that 
no appreciable interface dipole is formed. 

In absence of an interface dipole, the Schottky bar- 



rier at a MOI can be predicted by aligning the vac- 
uum levels of the metal and the organic material, called 
the Schottky-Mott rule. Over the last decade, however, 
experimental studies have indicated the general break- 
down of the vacuum alignment rule, and have demon- 
strate d that signifi cant interface dipoles are formed at 
^Qjgl8i22i2 3 | 24 | 25|26| addition such studies have shown 
that interface dipoles at MOIs are localized foremost at 
the first molecular layer covering the metal. The interface 
dipoles are not affected much by deposition of additional 
organic layers. Because MOI dipoles are localized at the 
interface, they can be extracted from the change in the 
surface work function after deposition of a single organic 
layer. 

Ideas inspired by chemical bonding have been put for- 
ward to explain large interface dipoles. Conventional 
semiconductors such as Si have reactive surfaces, which 
bind strongly to metal over layers. A significant density 
of states is then often created at the metal-semiconductor 
interface within the band gap of the semicon duct or, the 
so-called metal induced gap states (MIGS).^'^ '^MMi^ this 
model, MIGS determine the charge distribution at the in- 
terface and hence the interface dipole. The MIGS model 
has also been applied to MOIs.'^" '^^^^ It requires a strong 
interaction between the metal and the organic material. 

If molecules are physisorbed onto a metal surface, 
one expects a relatively weak interaction between the 
molecular semiconductor and the metal. For physisorbed 
molecules interface dipoles at MOIs have been explained 
by the so-called pillow effect .t^^"* ■'^'^ If a molecule 
approaches a metal surface, the electronic clouds of the 
molecule and the metal start to overlap. The Pauli ex- 
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change repulsion between these clouds leads to a spatial 
redistribution of electrons, which modifies the surface 
dipole. Since the electronic cloud of the metal is usu- 
ally "softer" than that of the molecule, the net effect of 
Pauli repulsion is that electrons are pushed back into the 
metal. The result is an interface dipole that decreases 
the work function of the surface. 

A decrease of the work function is commonly found if 
inert a toms or small m olecules are adsorbed on a metal 
surface .f^^lMSIinilSTlll] Remarkably, adsorption of larger, 
TT-conjugated, molecules can lead to a substantial in- 
crease, as well as a decrease of the work function and the 
dependence of this work function shift and the associated 
interface dipole on the molecules and the m etal has been 
the subject of intensive experimental stndyW^^^^^^^^^^^ 
If the work function W oi a, surface after coverage with 
a molecular layer is measured for a range of metal sub- 
strates with different initial work functions Wc, the re- 
sults can be characterized by the parameter 



where Wc, W are the work functions of the clean metal 
surface and of the surface with the adsorbed organic 
layer, respectively. 

The vacuum level alignment, or Schottky-Mott rule 
gives S = 1. Assuming that the pillow effect does not 
depend strongly on the metal substrate, it gives a rela- 
tively constant decrease of the work function, leading to 
5' « 1. Although this is observed for some molecules, 
very often S is significantly smaller than I.'* Moreover, 
there is no a priori reason why S should be a constant. 
Indeed for some molecules and polymers several regimes 
can be distinguished, between which a transition from 
5 w 1 to S* « is observed.tSaHa For the case where 
the organic layer is separated from the metal electrode 
by a thin insulating barrier, this behavior is interpreted 
with a model that assumes charge transfer across this 
barrier."*^ 

In this paper we study the dipole formation at in- 
terfaces of monolayers of PTCDA (3,4,9,10-perylene- 
tetra-carboxylic-di-anhydride), perylene and benzene 
molecules adsorbed on close-packed metal surfaces of Au, 
Ag, Al, Mg and Ca by means of density functional the- 
ory (DFT) calculations. We have selected these surfaces 
because they have a similar and simple structure, as well 
as a simple electronic (free electron like) structure. Yet 
their work functions span a range from 3.0 eV (Ca) to 5.3 
eV (Au), allowing to study the effect of the metal work 
function on the interface dipole. A preliminary account 
of this work has been given in Ref. Hil 

The molecules are chosen on account of their difference 
in complexity and their different behavior experimen- 
tally. PTCDA is a fairly complex conjugated molecule 
with a relatively small electronic gap. From an experi- 
mental point of view a PTCDA monolayer on metal sur- 
faces has been a model system to study MOIs. Deposi- 
tion of PTCDA on noble metal surfaces leads to well- 



ordered epitaxial over layers.!^ In particular the struc- 
ture and electronic stru cture of PTCDA on Ag (lll) have 
been studied intensively.Bl &7 | 48 | 49 | 5U | 5i | 5Li |53llllwork func- 
tion measurements have been performed for PTCDA ad- 
sorbed on a range of metal surfaces.* These mea- 
surements give a work function that is roughly indepen- 
dent of the metal substrate, i.e. S* « 0, meaning that ad- 
sorption of a PTCDA monolayer on a high work function 
metal gives a decrease of the work function, whereas ad- 
sorption on a low work function metal gives an increase. 

In contrast, experimental data suggest that adsorption 
of the simple conjugated molecule benzene on a metal 
surface always gives a decrease of the work function, with 
S in the range O.G-l.O.I^MIIIIIllIinil] ^he size and com- 
plexity of the perylene molecule is between that of ben- 
zene and PTCDA. The structure of an adsorbed perylene 
monola yer is thought to be similar to that of a PTCDA 
j^g^ygj. I48|62i63| Depending on the metal substrate, the work 
function can decrease or increase, but there is no uniform 
pinning as for PTCDA.ISH 

This paper is organized as follows. In the next sec- 
tion we give the technical details of our calculations. In 
Sec. |III| we present our results obtained for adsorption 
of PTCDA monolayers on different metal surfaces. We 
compare results obtained with different density function- 
al, and study the influence of the packing density of the 
molecules on the surface. Section |IV| gives the results 
obtained for adsorbed benzene and perylene monolayers. 
The results are discussed in Sec. |V] with the help of a 
simple phenomenological model, and a short summary 
and conclusions are given in Sec. |VI| 



II. COMPUTATIONAL DETAILS 

The electronic struct ure is treated within density 
functional theory (DFT^^EU ugjng the local density 
approximation (LDA),^^°^ or the generalized gradi- 
ent approximation (GGA) with the PW91 exchange- 
correlation functional.'^ The calculations are performed 
with the VAS P (Vienna ab initio simulation package) 
programJ ^° l ^^ l which uses the projector augmented wave 
(PAW) method. l22Esl For Au and Ag atoms the outer shell 
d and s electrons are treated as valence electrons, for Al 
the outer shell s and p electrons, and for Mg and Ca the 
outer shell s electrons. For atoms of first row elements 
the 2s and 2p electrons are treated as valence electrons. 
The valence pseudo wave functions are expanded in a ba- 
sis set consisting of plane waves. All plane waves up to a 
kinetic energy cutoff of 400 eV are included. 

To model the metal-molecule interface, we use a super- 
cell containing a slab of at least three metal layers with 
one layer of molecules adsorbed on one surface, and a 
vacuum region of at least 10 A. Periodic boundary con- 
ditions are applied and the atomic positions in the top 
metal layer and in the molecules are allowed to relax. A 
dipole correction is applied to avoid spurious interactions 
between dipoles of repeated slabs along the direction nor- 
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mal to the surface!^ 

The electronic structure is calculated self-consistently 
using a 3 X 3 (for PTCDA) to 5 x 5 (for perylene and 
benzene) k-point grid in the irreducible surface Brillouin 
zone (SBZ) according to the Monkhorst-Pack scheme 
and applying a Methfessel-Paxton smearing of 0.2 eV!^ 
A 3 X 3 k-point grid gives well-converged results for 
PTCDA layers, because of the large size of the surface 
unit cell (see next section). For accurate calculations of 
total energies and densities of state (DOS) the charge 
densities are recalculated with an up 7 x 7 k-point grid, 
using the tetrahedron method.l^^ DOSs are plotted using 
Gaussian smearing with a broadening parameter of 0.1 
eV. 

Work functions are evaluated from the expression 

W ^Vioo) - Ef, (2) 

where V{oo) is the electrostatic potential in the vacuum 
region and Ep is the Fermi energy of the bulk metal. 
V{oo) is obtained from the potential averaged in the 
(x, y) plane 

^(^) = 4 // V{x,y,z)dxdy, (3) 

^ J J cell 

where V{x, y, z) is the electrostatic potential on a real 
space grid in the supercell. In practice V{z) reaches an 
asympt otic v alue V{(X)) at a distance of a few A from the 
surface .^2*211 An accurate value of Ep is obtained from 
a separate bulk calculation, following the procedure de- 
scribed in Ref. ,78, 

Test calculations regarding slab thickness, vacuum 
thickness, k-point sampling grid and plane waves kinetic 
energy cutoff are performed in order to estimate the con- 
vergence. From these tests we find that with the param- 
eters given above, total energies are converged to within 
0.01 eV and work functions to within 0.05 eV. The results 
for PTCDA on Ca(lll) turn out to be the most sensitive 
with respect to vacuum thickness and k-point sampling. 
So for this system the results have been obtained using a 
5x5 k-point grid in the irreducible SBZ and a vacuum 
thickness of at least 14 A. 

To analyze our results we also use properties of isolated 
molecules, such as the electron affinity (EA), as will be 
discussed in Sec.|Vl For calculations on isolated molecules 
we use the GAMESS program,E2l and treat the elect ronic 
structure within DFT using the BLYP functionaLl^SlSil 
We use the 6-31+G* basis set, which gives EAs for acenes 
that are converged on a scale of ~ 0.1 eV.I^ As in Ref. [5^ 
we find that including a diffuse orbital in the basis set is 
important and that the smaller 6-31G* basis set does 
not give a sufliciently converged EA.^'^ For instance, the 
EAs of PTCDA obtained using the 6-31G* and 6-31-HG* 
bases differ by 0.4 eV. The Kohn-Sham energy levels of 
the (neutral) isolated molecule calculated with VASP and 
GAMESS are very similar, illustrating that, provided the 
basis sets are sufficiently converged, the PW91 and BLYP 
functionals give similar results. 



TABLE I: Optimized nearest neighbor distances in the bulk 
metals. All values are in A. 





Au 


Ag 


Al 


Mg 


Ca 


GGA 


2.94 


2.93 


2.86 


3.19 


3.92 


LDA 


2.87 


2.84 


2.82 


3.13 


3.78 


exp. 


2.88 


2.89 


2.86 


3.21 


3.95 



TABLE II: Calculated work functions of clean (111) surfaces; 
(0001) for Mg . All values are in eV. 





Au 


Ag 


Al 


Mg 


Ca 


GGA 


5.25 


4.50 


4.08 


3.74 


2.98 


LDA 


5.52 


4.90 


4.21 


3.93 


3.08 


exp. 


5.26," 5.35* 4.46," 


4.50,'' 


4.56" 4.24^ 


3.78" 


(2.87)' 


calc. 


5.27,' 5.35* 


4.42^ 


4.25* 


3.76,' 3.8: 


2.86" 



"Ref. mi 
''Ref. [13 
'^Ref. |85] 
''Ref. in 
"Ref. 113 
f'Rei. [88] 
fRef. |89l 

''polycrystalline value, Ref. 1901 
'GGA, Ref. El] 
JGGA, Ref.|92] 
'■XDA, Ref. [93] 
'LDA, Ref.l94l 
"GGA, Ref.[9il 
"LDA, Ref. ESI 

III. PTCDA 

Before discussing the results obtained for adsorbed lay- 
ers, we benchmark our calculations on clean metal sur- 
faces. We consider the close-packed (111) surfaces of fee 
Au, Ag, Al, and Ca, and the (0001) surface of hep Mg. 
The metals in this set are relatively simple, free electron 
like and the set of surfaces spans a considerable range 
in work functions. Table |l] lists the optimized nearest 
neighbor distances of the bulk metals, calculated with 
GGA(PW91) and LDA functionals. As usual, the GGA 
values are larger than the LDA values, but both are gen- 
erally in reasonable agreement with experiment. We use 
these optimized distances to construct the surface unit 
cells. 

Table |TT] lists the calculated work functions of the clean 
(111) surfaces (for Mg the (0001) surface), compared to 
experimental values and values obtained in previous cal- 
culations. Our results have been obtained using slabs 
consisting of six metal layers. A 25 x 25 k-point sampling 
of the SBZ is applied, while allowing the top two metal 
layers to relax. The GGA values generally are within 
^ 0.1 eV of the experimental values. LDA still gives an 
acceptable accuracy, but tends to overestimate the work 
function somewhat. Only for Al(lll) LDA gives a bet- 
ter value than GGA, as compared to experiment. Our 
results also agree with those obtained in previous com- 
putational studies; the small differences can be attributed 
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FIG. 1: (Color online) PTCDA monolayers on the Ag(lll) 
surface. The rectangles denote the surface unit cells used in 
calculations; (a) the dilute structure (with area A = 268A^); 
(b) the herringbone structure {A — 243A'^). 



to differences in the computational parameters, such as 
the functional, the basis set, and the lattice parameter. 



A. Structure of adsorbed monolayers 

PTCDA monolayers adsorbed on Ag(lll) and 
Au(lll) surf aces have been studi ed in detail 
experimentally! 4" l ^7 | 4'^ l ^" l ^" l ^^ l ^" l ^^ l "'^ l »7 | "^ l In a close- 
packed monolayer the PTCDA molecules lie flat on the 
surface in a "herringbone" structure with the centers 
of the PTCDA molecules located on surface bridge 
sites.!^ The surface unit cell contains two PTCDA 
molecules, see Fig. [ijb). The experimental distances 
between the carbon rings of the molecules and the 
surface atoms are 2.86 A and 3.27 A for ad sorption 
on Ag(lll) and Au(lll), respectively.!^^!^ Experi- 
ments indicate a weak interaction betwe en P TCDA 
and Au(lll), consistent with physisorption,^^'^!] and a 
somewh at stronger interaction between PTCDA and 
Ag(lll).Elll |4»|5U | 5i | PTCDA binds more strongly to 
open Ag surfaces and to surface steps .1^22! We do not 



TABLE III: Average bond lengths (in A) and bond angle of 
PTCDA adsorbed on Ag(lll). 





LDA 


CCA 


GGA 




dilute 


dilute 


herringbone 


C-H 


1.10 


1.09 


1.09 


c-c 


1.42 


1.43 


1.43 


C-O" 


1.22 


1.23 


1.23 


C-0'' 


1.38 


1.40 


1.40 


-0-C(°) 


125.2 


125.3 


124.7 



"carboxyl 
'anhydride 



know of any such detailed studies on the structure of 
PTCDA adsorbed on the other metal (111) surfaces. 
We refrain from comparing our results to experiments 
where metals are deposited onto thin films of PTCDA, 
as this often lead s to interdiffusio n, which complicates 
the interpretation! i"i | i"^ | i"^ | i"4 | i 05l 

In our calculations we use the herringbone structure of 
PTCDA on Ag shown in Fig. [ijb). The underlying Ag 
substrate contains 33 metal atoms per unit cell and our 
supercell contains 175 atoms in total. Since the lattice 
parameters of Au, Ag and Al are similar, see Table |T] 
we use a similar supercell for PTCDA on these surfaces. 
For Mg and Ca, we choose a herringbone structure that 
results in a packing density of PTCDA molecules similar 
to that on the other surfaces. This results in 30 and 20 
atoms per metal layer for Mg and Ca, respectively. 

To study the effect of the packing density of PTCDA 
molecules, we also perf orm calculations on the structure 
used by Picozzi et al.^^ see Fig. [ija). This surface 
unit cell contains one PTCDA molecule. We refer to 
this structures as the "dilute" structure. The under- 
lying Ag(lll) surface then contains 36 atoms per unit 
cell, so that the coverage of PTCDA molecules is ^ ^ 
ML. The surface unit cells for the other metal substrates 
are chosen such, that the coverage remains close to this 
value. The distance between the PTCDA molecules is 
then sufficiently large for the molecules to have no di- 
rect interaction. As it is easier to vary the geometry 
of the molecule and substrate in the dilute structure, as 
compared to the close-packed herringbone structure, we 
use the former to study the energetics of PTCDA adsorp- 
tion. The optimized geometries of the PTCDA molecules 
in the two structures are very similar, as demonstrated 
by Table |III| suggesting that close packing the molecules 
in the herringbone structure does not lead to a large in- 
termolccular interaction. In addition, the GGA or LDA 
optimized geometries are very similar. 

Common DFT functionals describe strong chemical 
interactions well, but they fail to capture weaker (van 
der Waals) bonding correctly. If GGA functionals are 
used to describe the adsorption of closed-shell molecules 
on metal surfaces, this can lead to underestimating the 
binding energy and overestimating the bond distance 
between molecule and surface. Using LDA functionals 
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FIG. 2: (Color online) Binding energy curves for (planar) 
PTCDA on Au(lll) (squares) and Al(lll) (circles). GGA 
and LDA values are indicated by the dashed and solid lines, 
respectively. 



however can lead to a serious overbmding and a equilib- 
rium distance that is too small'^^'^"^' Previous GGA cal- 
culations of the binding energy of PTCDA on Ag(lll) 
give slightly varying results, i.e. a moderate binding 
of ^ 0.5 eV/molecule,'^^ or a very weak binding of 
< 0.1 eV/ molec ule, or even a purely repulsive binding 
curve) 

WM^imim In the calculations where binding 
was obtained, the equilibrium distance ('^ 3.4 A) is signif- 
icantly larger than the experimental equilibrium distance 
(2.9 A). 1^21531 contrast, LDA calculations on PTCDA 
on Ag(lll) give equilibrium distances of 2.8 A^°^ and 
2.7 A P which are in better agreement with experiment. 
The LDA binding energy, ~ 3 eV/molecule, however, is 
suspiciously large.l^ We obtain very similar results in our 
GGA and LDA calculations for PTCDA on Ag(lll). 

To illustrate how general this trend is. Fig. [2] shows 
binding energy curves calculated with the dilute struc- 
ture , where we varied the distance between the PTCDA 
molecule and the Au and Al(lll) surfaces. The GGA 
results for PTCDA on Au lead to an extremely shallow 
binding curve with a minimum at a distance > 4 A and 
a very small binding energy of ~ 0.1 eV/molecule. Using 
GGA for PTCDA on Al gives a sizable binding energy 
('~ 0.8 eV/molecule) and an equilibrium distance ~ 3.5 
A. LDA calculations lead to much larger binding ener- 
gies, i.e. 1.7 eV/molecule for PTCDA on Au and 3.0 
eV/molecule for PTCDA on Al. The corresponding equi- 
librium distances are 3.15 A and 2.95 A, respectively. An 
LDA calculation for PTCDA on Ag gives an equilibrium 
distance of 2.75 A, which is somewhat smaller that the 
experimental value. 

In conclusion, GGA and LDA give different results 
for the binding in weakly bonded systems. GGA gives 
small molecule-surface binding energies and large equi- 
librium distances, whereas LDA gives large binding ener- 
gies and smaller equilibrium distances. As the reactivity 



FIG. 3: (Color online) Optimized geometry of the arching 
PTCDA molecule on the Ca(lll) surface. 



of the surface increases along the series Au, Ag, Al, Mg 
and Ca, one expects the adsorption of PTCDA gradually 
changing from physisorption to chemisorption. The dif- 
ferences between the GGA and LDA become smaller and 
the results more rehable. For PTCDA on Mg(OOOl) and 
Ca(lll) we obtain GGA binding energies of 2.3 and 8.4 
eV/molecule, respectively, indicative of chemical bond- 
ing. 

We also find that, if the binding energy increases, the 
PTCDA molecule looses its planar geometry. If the in- 
teraction between the PTCDA molecule and the surface 
is large, the molecule arches as shown in Fig. |3] The 
extend of this geometry deformation depends upon the 
metal substrate. Optimizing the geometry of PTCDA 
on Ca with GGA and defining the position of the sur- 
face by the average z coordinate of the top layer of Ca 
atoms, we find that the outer carbon atoms of the pery- 
lene core are 0.8 A closer to the surface than the carbon 
atoms in the center. The latter have a distance of 2.6 A 
to the surface. Such short distances are indicative of a 
strong interaction between molecule and surface. A ge- 
ometry deformation also occurs in the end groups of the 
PTCDA molecule, where the carboxyl oxygens move to- 
wards the surface, and the anhydride oxygen moves way 
from the surface. A strong deformation of the PTCDA 
molecule is accompanied by a rumpling of the surface, 
where metal atoms are lifted out of the surface, decreas- 
ing the distance with the molecule. For instance, the 
distances between the carboxyl oxygens and the nearest 
Ca atoms are 2.3 A. Such short distances suggest the for- 
mation of bonds, which matches the large binding energy 
of 8.4 eV/molecule. 

The deformation of the adsorbed PTCDA molecule 
and that of the metal substrate decrease through the se- 
ries Ca, Mg, Al, Ag, and Au, and is accompanied by a 
decrease of the binding energy. The deformation pattern 
of the PTCDA molecule qualitatively remains the same 
along this series, but the amplitude of the deformation 
decreases. For instance, the (GGA) distances between 
the carboxyl oxygens and the nearest Mg atoms are 2.4 
A and the outer carbon and core carbon atoms are at 2.6 
and 2.8 A from the surface, respectively. This matches 
the binding energy of 2.3 eV/molecule, which is much less 
than between PTCDA and Ca(lll). Throughout the se- 
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ries Al, Ag and Au the binding energy, as well as the 
geometric deformation of the molecule, decrease mono- 
tonically. 

For instance, if we fix the overall PTCDA-Ag(lll) dis- 
tance at 2.75 A and optimize the geometry, the molecule 
is only slightly arched. The outer carbon atoms are 0.1 
A closer to the surface than the core carbon atoms. The 
carboxyl oxygen atoms are 0.2 A closer to the surface, 
whereas the anhydride oxygens are at approximately the 
same height as the core carbon atoms. This pattern is 
in fair agreement with experimental observations.-- At 
the end of the series, i.e. for PTCDA on Au(lll), the 
binding energy is vanishingly small and the molecule and 
surface have an undistorted, planar geometry. The bind- 
ing energies and geometries of PTCDA on Al and Ag 
obtained in previous calc ulations follow the trends dis- 

cussed above.fsSEnmniiinii 



B. Work functions 

From the calculations discussed in the previous sec- 
tion we can understand the trends in the bonding and 
in the geometry of PTCDA adsorbed on the different 
metal surfaces. In this section we study the work func- 
tion changes induced by the adsorption, which originate 
from a redistribution of charge at the MOI. The charge 
distribution is mainly determined by electrostatic and 
short-range exchange correlation interactions, which are 
represented well by LDA and GGA functionals) ^^° I ^^ ^For 
example, the charge transfer between the weakly inter- 
acting molecules in a charge transfer crystal is described 
well by LDA,'ii2J as is th e work function graphene ad- 
sorbed on a metal surface .'^^^'^^ Given the bonding ge- 
ometry of the molecule and surface, we expect therefore 
these functionals to give good values for the work func- 
tion of molecular monolayers adsorbed on a metal sur- 
face. However, since they do not incorporate van der 
Waals interactions correctly, calculations lead to a con- 
siderable uncertainty in the molecule-surface equilibrium 
distance for weakly bonded systems, as discussed in the 
previous section. Because of this uncertainty we inves- 
tigate the formation of interface dipoles in a number of 
steps. We start with the dilute structure and perform 
calculations for fixed molecule surface distances d in the 
range 3.0-3.6 A. Results obtained with GGA and with 
LDA are then compared. In the second step we switch 
to the more densely packed herringbone structure that 
is observed experimentally, which allows us to study the 
effect of the packing density. Finally, we discuss the ef- 
fects of full geometry relaxation of the molecules and the 
surface. 

Fig. |4]ja) shows the work functions for a layer of pla- 
nar PTCDA molecules adsorbed in the dilute structure 
on the different metal surfaces, calculated using the GGA 
functional. One immediate observation is that adsorption 
on Au(lll) leads to a lowering of the work function as 
compared to the clean surface, whereas adsorption on the 
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FIG. 4: (Color online) Work functions of a PTCDA monolayer 
on (111) metal surfaces [(0001) for Mg] in the dilute structure, 
calculated using (a) the GGA and (b) the LDA functionals. 
The clean metal work function is given along the i-axis. The 
numbers give the calculated values, the lines guide the eye. 
The bottom (black) curve refers to the clean metal surfaces. 
The top three curves are for different distances d between the 
PTCDA molecules and the surfaces. 



other metal surfaces leads to an increase of the work func- 
tion. There is some dependence of the work function on 
the distance between the molecule and the surface, but it 
is not excessively large. By fitting a straight line through 
the curves in Fig. [4]ja) one obtains (see Eq. ([T])) S* = 0.5 
at c? = 3.6 A, and S = 0.6 at d = 3.0 A. These values 
are considerably lower than the S = 1 that follows from 
the Schottky-Mott rule, indicating that significant inter- 
face dipoles are formed upon adsorption. Since the work 
function changes decrease somewhat upon decreasing the 
molecule-surface distance, the interface dipoles decrease 
with decreasing distance. The calculated values for S are 
much higher than the S ~ Q obtained experimentalljl^'^SI^ 
We will show below that this discrepancy is resolved by 
increasing the packing density of the PTCDA molecules, 
which is only ^ ^ ML in the dilute structure. 

To illustrate the effect of using a different functional, 
Fig. |4|b) gives the work functions for the dilute struc- 
ture at fixed distances between the molecule and the sub- 
strates, calculated using the LDA functional. Compared 
to the GGA results of Fig. Qa), the LDA work func- 
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TABLE IV: LDA optimized molecule-surface distances dmin 
and GGA work functions W of PTCDA monolayers in the her- 
ringbone structure, adsorbed on metal(lll) (for Mg (0001)) 
surfaces. 







W eV 


Au 


3.15 


4.78 


Ag 


2.75 


4.65 


Al 


2.95 


4.68 


Mg 


2.7-2.8 


3.80 


Ca 


2.3-2.4 


3.39 



tions are generally somewhat higher, as was also the case 
for the clean metal surfaces, see Table [ill The changes in 
the work functions upon adsorption calculated with LDA 
or GGA, are comparable however at the same molecule- 
surface distance. It means that, although the binding 
between the molecules and the surfaces calculated with 
LDA or GGA can be considerably different, as discussed 
in the previous section, the charge redistribution upon 
adsorption is similar, if we consider the same molecule- 
surface distance. Since the GGA work functions of the 
clean metal surfaces are somewhat closer to the experi- 
mental values, see Table [IT] we will use GGA values for 
the adsorbed layers in the following. 

Fig. [5] shows the work functions of the herringbone 
structure of PTCDA on metal surfaces, calculated us- 
ing the GGA functional. As for the dilute structure. 
Fig. Qa), adsorption on Au(lll) lowers the work func- 
tion, and on the other metal surfaces it increases the 
work function. In the herringbone structure the work 
function shifts are much larger however. Most strikingly, 
as one can observe in Fig. [S] the work function is pinned 
at ~ 4.7 eV over a considerable range of metal substrates 
and molecule-surface distances. The pinning leads to 
S = 0, which is in agreement with experiment. More- 
over, the value of the pinned wor k function is close to 
the value found experimentallj^SEH^ Deviations from pin- 
ning are observed only for low work function metals and 
short molecule-surface distances, i.e. d < 3.3 A for Ca 
and d = 3.0 A for Mg. 

LDA gives very short molecule-surface distances, which 
can serve as lower bounds (i„iin- We have optimized the 
geometries at dmin of an adsorbed PTCDA monolayer 
in the herringbone structure and calculated the work 
functions with GGA, see Table |IV| These work func- 
tion results for adsorption on Au, Ag and Al(lll) are 
very close to the pinning value in Fig. [5] The work 
functions for adsorption on Mg(OOOl) and Ca(lll) drop 
to a lower value. As discussed in the previous section, 
the Ca and Mg surfaces become somewhat unstable at 
dmin, i-e. metal atoms are pulled up from the surface, 
which is accompanied by a strong arching distortion of 
the PTCDA molecule. Such a deformation generates a 
molecular dipole moment perpendicular to the surface 
that decreases the work function. 
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FIG. 5: (Color online) As Fig. [i] for PTCDA in the herring- 
bone structure 



IV. BENZENE AND PERYLENE 

Although there is some spread in the numbers mea- 
sured, experiments clearly indicate a lowering of the 
work function after adsorption of a benzene monolayer 
on several metal surfaces. Measured work function 
shifts ar e 0.18 e V on A1(111),E5I 0.3 eV and 7 eV on 
Ag(lll)IIZIlllI]o.7 eV and 1.05 eV on Cu(lll),l23l2land 
1.10 eV on Au(lll).'23 Previous GGA calculations for 
benzene on Al(lll) gave an equilibrium distance of 3.7- 
3.8 A.E2I Quantum chemical MP2 calculations for a sin- 
gle benzene molecule adsorbed on a cluster of metal 
atoms gave equilibrium distances of 3.8 A for benzene on 
Au(lll) and 4.0 A for benzene on Cu(lll)E3. We calcu- 
late the work functions of a benzene monolayer adsorbed 
on different metal surfaces at a set of fixed distances in 
the same way as discussed in the previous section. A 
(\/7 X \/7)-R19.1° structure is used as in Ref. [B51 which 
gives a somewhat less than close-packed coverage of the 
surfaces. 

The (GGA) results are given in Fig. [e] The most im- 
portant observation is that adsorption of benzene leads 
to a decrease of the work function for all the surfaces 
studied. This is very different from the effect of PTCDA 
adsorption, see Figs. |4] and [5] Moreover, adsorption of 
benzene gives a work function lowering that is of a sim- 
ilar size for all surfaces (at a fixed molecule-surface dis- 
tance). This leads to S" = 0.9 at d = 3.6 A, and S = 0.8 at 
d — 3.0 A, which is not extremely far from the Schottky- 
Mott limit S = 1. The absolute size of the work function 
shift depends on the molecule-surface distance with the 
number at d = 3.0 A being roughly twice as large as 
that at d = 3.6 A. The sign of the work function shift, 
its relatively weak dependence on the metal, and its sen- 
sitivity to the molecule-metal distance all point to an 
interpretation in terms of the pillow effect. The effect is 
determined by the Pauli repulsion between the molecular 
and surface electrons, which dec reases th e surface dipole 
and therefore the work function .'22I2MI11 Pauli repulsion 
critically depends on the overlap between the molecular 
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FIG. 6: (Color online) As Fig. [4] for benzene in a (\/7 
V7)R19-1° structure.-^ 



and surface wave functions and therefore on the distance 
between the molecule an d the surface. Our calculations 
and previous calculation^'^^'^^' suggest that the distances 
between the benzene molecule and the metal surfaces are 
rather large, i.e. > 3.5 A. 

Summarizing, adsorption of the large PTCDA 
molecule leads to work function pinning (S ~ 0), and 
adsorption of the small benzene molecule gives a uniform 
work function lowering (S' « 1). It is then interesting to 
study the adsorption of an intermediately sized molecule, 
such as perylene. The structure of a perylene monolayer 
on a metal surface is less well-established than that of 
a PTCDA monolayer. A herringbone structure similar 
to PTCDA is the structures proposed for a close-packed 
perylene monolayer on Ag(lll) and Au(lll).^® ''■^ UPS 
measurements give a decrease of the work function of 
perylene on Au and Ag by 0.8 eV and 0.6 eV, respec- 
tively, and an increase of the work function of perylene 
on Ca by 0.3 eV.EH 

We perform calculations for a perylene monolayer ad- 
sorbed on different metal surfaces at a set of fixed dis- 
tances. Analogous to PTCDA we use two different struc- 
tures, i.e. a close-packed herringbone structure, and a di- 
lute structure with a packing density of ML. The cal- 
culated work functions are shown in Fig. u\ Two regimes 
can be distinguished. For the high work function surfaces 
(Au, Ag), adsorption of perylene leads to a lowering of 
the work function, whereas for the low work function 
surface of Ca, adsorption of perylene increases the work 
function. These results are in qualitative agreement with 
experiment. The curves in Fig. [7] indicate that the tran- 
sition between these two regimes takes place in the range 
Mg-Al. 

From these curves the transition between the two 
regimes can be quantified. Starting with the results ob- 
tained for the dilute structure at a molecule-surface dis- 
tance d = 3.6 A, a line through the points for Ca, Mg 
and Al gives S = 0.3, whereas a line through the points 
for Al, Ag and Au leads to 5" = 1.0. For the herring- 
bone structure the same procedure for d = 3.6 A gives 
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FIG. 7; (Color online) As Fig.[4]for perylene in (a) the dilute 
structure and (b) the herringbone structure. 



5 = and 5* = 0.9, respectively. It is instructive to 
compare these S values to the values obtained for ben- 
zene and PTCDA. It suggests that for Ca, Mg and Al 
one obtains pinning of the work function upon perylene 
absorption, similar to PTCDA, see Fig. [s] whereas for 
Ag and Au one finds a uniform work function decrease, 
similar to benzene, see Fig. [6] Upon decreasing the dis- 
tance between the perylene molecules and the surfaces 
the S values in the high and low work function regimes 
become somewhat closer and the transition between the 
two regimes becomes less sharp. Note that the distance 
dependence in the low work function regime resembles 
the distance dependence of the PTCDA case, whereas in 
the high work function regime it resembles the benzene 
case. Experimentally such a transition is observed for 
Alqa adsorbed on different surfaces. For adsorption on 
low work function metals 5 « 0, whereas for adsorption 
on high work function metals 5 « l.'22l 



V. DISCUSSION 

The results presented in Fig. |5] show that adsorption 
of a PTCDA monolayer pins the work functions at ~ 4.7 
eV for a broad range of distances. The fact that one gets 
pinning, as well as the value of the pinning level, are in 
good agreement with experimental observations. The re- 
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suits obtained for benzene and perylene adsorption are 
in qualitative agreement with available experimental re- 
sults, i.e. adsorption of benzene leads to a lowering of 
the work function in all cases, and adsorption of pery- 
lene gives a work function decrease for high work func- 
tion metals and a work function increase for low work 
function metals. In this section we analyze this behavior 
and interpret the results. First we analyze the density of 
states, then we consider explicitly the charge transfer at 
the molecule-substrate interface, and finally we formulate 
a simple phcnomcnological model. 



A. Density of states 

Figure [s] gives the Kohn-Sham (KS) density of states 
(DOS) of an isolated PTCDA molecule, calculated us- 
ing the GGA functional. The energies of the highest 
occupied molecular orbital (HOMO) and the lowest un- 
occupied molecular orbital (LUMO) are at —6.24 eV re- 
spectively —4.64 eV with respect to the vacuum level. 
The KS spectrum is similar to that obtained in previous 
DFT calculations.l^nEniiniinilThe GGA HOMO-LUMO 
gap of 1.60 eV also agrees with the value found in other 
GGA calculations.30 

The interpretation of KS energy levels is the subject of 
a long-standing debate. From calculations with advanced 
functionals it is argued that the KS energies of all occu- 
pied molecular orbitals correspond to vertical ionization 
potentials (IPs), as can be extracted from a photoemis- 
sion spectrum, for instance. Approximative functionals 
such as GGA then still give a reasonable ionization spec- 
trum, but it is s hifted t o a higher energy by approxi- 
mately a const ant .l^iSinnixhe work function is the lowest 
IP of an extended system. Even approximative function- 
als such as GGA or LDA usually give work functions 
that are close to the experimental values, as is illus- 
trated by Table |ll] Results of a similar quality are ob- 
tained for work functions of adsorbed atomic and molec- 
ular layers.EsEEIIIMnili] 

The KS energy levels corresponding to unoccupied 
molecular orbitals generally do not have such a simple 
interpretation. In particular, the energy of the DFT 
LUMO (eg) should not correspond to the electron affin- 
ity (EA). From calculations with accurate functionals it 
is shown that eq < — E A, bo th for molccules,'^^^! as well as 
for extended systems! "*^^^*^^^ The difference between |eo| 
and the EA can be several eV's, which is attributed to 
the fact that an accurate functional has a discontinuous 
derivative as function of the number of electrons. 

A similar difference is found for approximative con- 
tinuous functionals such as GGA or LDA. However, it 
is well-known that for these approximative functionals. 
Slater's transition state approach allows for a simple es- 
timate of the E A .1^221 "We define eo as the KS energy of 
the LUMO of the neutral molecule, and eM(l) as the KS 
energy of the singly occupied HOMO of the ion that has 
one additional electron. The functionals allow for a frac- 



tional occupancy of t he le vel ej\/ with N electrons, and 
using Janak's theorenP^l 



dN ' 



where i?tot is the total energy, one can write 



EA 



£M{N)dN. 



(4) 



(5) 



To a good approximation the dependence of on the 
occupancy is linear and can be parameterized as 



(6) 



where U is the effective charging energy per 
electron .E^IElll This then gives 



EA = -eo - ^U. 



(7) 



This procedure gives results that agree very well with 
charging energies for isolated conjugated molecules ex- 
tracted from total energy calculations.^*' 

From separate SCF calculations on the neutral 
molecule and the singly charged ion one can extract eo 
and eM(l), and calculate U = eM(l) — eo- Following this 
procedure, we extract a charging energy U = 3.31 eV 
from calculations on PTCDA° and PTCDA". Using the 
LUMO energy eq = —4.64 eV, we then find from Eq. ([T]) 
EA = 2.98 eV for PTCDA. This value is in very good 
agreement with the value EA = 2.96 eV we extract from 
a ASCF total energy difference calculation. Slater's tran- 
sition state approach can also be used to calculate the IP. 
Assuming that the charging energy for holes on PTCDA 
is the same as for electrons and using the HOMO energy 
of —6.24 eV, then yields IP = 7.90 eV, which is in fair 
agreement with the experimental PTCDA gas phase IP 
of 8.15 eV cited in Ref. [TTCI 

If a molecule is embedded in a crystal, its charging 
energy U is reduced drastically because of screening by 
the surrounding crystaL'^^'^^Zl One can write 



C/= C/barc-2P-, 



(8) 



where Uhaie is the charging energy of a bare molecule, and 
P~ is the polarization energy associated with a singly 
charged molecular ion in a crystal.l^^Il Using the polariza- 
tion energy P~ — 0.91 eV from Ref. 11271 we obtain an 
effective charging energy of a PTCDA molecule embed- 
ded in a PTCDA crystal U„yst = 3.31 - 1.82 = 1.49 eV. 



cryst 

Slater's transition state model then gives EA, 



cryst 



3.90 



eV and IPcryst = 6.99 eV for the EA and IP of a PTCDA 
molecule in the crystal. This leads to a transport gap 
Et ^ IPcryst - EAcryst = 3.09 eV. The values of the 
IP and Et are in good agreement with the values ex- 
tracted from ex perime nt, i.e 6.7 ± 0.2 eV and 3.2 ± 0.4 
eV, respectively !^^ ' ^^^! 

In summary, using the KS DOS as shown in Fig. |8] 
to calculate measurable quantities, one has to incorpo- 
rate the charging energy U of the molecule, cf. Eqs. ^ 
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FIG. 8: Total DOS of an isolated PTCDA molecule calcu- 
lated using a Gaussian broadening of 0.01 eV. The highest 
occupied molecular orbital (HOMO) and the lowest unoccu- 
pied molecular orbital (LUMO) are indicated. The energy is 
set to zero at the position of the HOMO. With respect to the 
vacuum level the HOMO and LUMO levels are at -6.24 eV 
and —4.64 eV, respectively. 



and ^ . U strongly depends upon the interaction of the 
molecule with its environment. Screening by a me tal sub- 
strate reduces U significantly, for instance P^^ * ^^^ * and for 
electrons in completely delocalized states U — 0. 

We now consider the DOS of PTCDA monolayers ad- 
sorbed on metal surfaces. Since the full DOS is domi- 
nated by states originating from the metal substrate, we 
look at the DOS projected on the atoms of the molecule 
(PDOS), in order to identify the contributions of the 
molecular levels. Fig. |9] gives the PDOSs PTCDA mono- 
layers adsorbed on metal surfaces at distances of d — 3.0 
and d = 3.6 A. We have applied a Gaussian broadening 
with a broadening parameter of 0.1 eV to avoid a spiky 
appearance of the PDOS. Comparison to Fig. |8] allows 
us to identify the the contributions of the molecular lev- 
els, in particular those states that are dominated by the 
molecular HOMO and LUMO. Although all states result 
from hybridization between molecular and metal states, 
it simplifies the analysis if we label them by their domi- 
nant molecular character. 

The interaction between the molecule and the sur- 
faces induces hybridization of molecular and metal states. 
For adsorption on simple, wide band metals this leads 
to a simp le b roadening of the molecular levels into 
resonances .'i^ For adsorbed PTCDA, the extend of this 
broadening is moderate. At d = 3.6 A the typical width 
at half height of the HOMO and LUMO peaks is ~ 0.2 
eV. The widths increase with decreasing molecule-metal 
distances (to ~ 0.3 eV at d = 3.0 A, for instance). These 
results agree with the widths found in previous DFT cal- 
culations of adsorbed PTCDA, as well as with those 
typically foun d for other adsorbed molecules such as 
pentacene.'^^They are also in qualitative agreement with 
peak widths observed in scanning tunneling spect roscopy 
of PTCDA adsorbed on Ag and Au(lll) surf aces. l^iEMS] 
The widths found in the calculations of Refs. [301 and |3T] 




FIG. 9: Projected density of states (PDOS) of the PTCDA 
molecule adsorbed on metal surfaces at a fixed distance d in 
the herringbone structure, calculated using a Gaussian broad- 
ening of 0.1 eV; left: d = 3.0A; right: d = 3.6A. The peaks 
corresponding to the molecular HOMO and LUMO levels are 
labeled. 



are much larger. 

Comparing the two columns in Fig. [9j one observes 
that the PDOSs at the two distances are qualitatively 
similar, but that the spectrum at d = 3.0 A is shifted 
towards lower energy as compared to the spectrum at 
d — 3.6 A by up to 0.5 eV, depending upon the metal 
substrate. This is caused by the pillow effect, as we will 
discuss in Sec. |VC] The GGA HOMO-LUMO gap of 
PTCDA decreases somewhat upon adsorption, from 1.6 
eV in the isolated molecule to 1.4 eV (peak maximum 
to peak maximum) for PTCDA weakly interacting with 
Au(lll) and ~ 1.1 eV for PTCDA strongly interacting 
with Ca(lll). 

One observation that can be made by comparing 
Figs, [s] and |9] is that work function pinning occurs 
when the Fermi level crosses the level of the LUMO. For 
PTCDA on Au the LUMO is unoccupied, but already 
for PTCDA on Ag the LUMO gets partially occupied. 
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This implies that electron transfer takes place from the 
metal substrate to the molecule. Judging from the up- 
wards shift of the Fermi level along the columns of Fig. [9] 
the amount of electron transfer increases along the se- 
ries Ag, Al, Mg, and Ca. As long as the Fermi level 
is inside the LUMO peak, the work function is pinned, 
compare Fig. [5] At the short molecule-surface distance 
of 3.0 A between PTCDA and Ca, the Fermi level jumps 
to the next peak, i.e. the LUMO-t-1, see Fig. [oji). This 
is accompanied by an "unpinning" of the work function, 
compare Fig. [5] 

The most detailed experiments have been performed 
for PTCDA on Ag(lll). In UPS and STM experiments 
a peak is observed at the Fermi level that is identified 
as the LUMO of the PTCDA molecule, whereas a peak 
at ~ —1.6 eV with resp ect to the Fermi level is labeled 
as the HOMO .1^21511521551 pjg jgjj^-)^ .jyijgj.g ^Yie distance 

between PTCDA and the Ag surface is close to the ex- 
perimental value, we find the HOMO at —1.3 eV and the 
LUMO at the Fermi level, which is in reasonable agree- 
ment with the experimental analysis. In STM and IPES 
experiments of PTCDA on Au(lll) a peak is observed at 
1.0-1.5 eV above the Fermi level that is associated with 
the EA level of the molecule .12111231111. Applying Eqs. Q 
and Jsl puts the EA level at P~ above the LUMO in 
Fig.p[a,b). Using the polarization energy P~ = 0.97 eV 
for a PTCDA ion adsorbed on a metal surface calculated 
in Ref. I127l then brings the EA level in the experimentally 
observed range. 

Pinning of the work function at MOIs has bee n inter- 
preted in terms of a charge neutrality level (CNL) po | 3i | 32 | 
in analogy to Schottky b arrier models for conventional 
semiconductors .I^SESHIMI xhe CNL model relies on hav- 
ing a large continuum DOS at the metal-semiconductor 
interface, which fills the energy gap of the semiconductor. 
The Fermi level is then pinned by these metal induced 
gap states (MIGS).'23 Conventional semiconductors such 
as Si or GaAs have reactive surfaces with surface atoms 
carrying dangling bonds. The energies of these dangling 
bond states are within the semiconductor gap. Bonding 
at a metal-semiconductor interface leads to broadening 
of these states, which generates a large continuum DOS 
at the interface in the semiconductor ga,p^^^^^ Closed 
shell molecules such as PTCDA do not have dangling 
bond states within the HOMO-LUMO gap. The creation 
of a large DOS at a MOI within the HOMO-LUMO gap 
then depends upon a large broadening of the molecular 
levels. We do not observe such a large broadening. 

The DOS of an isolated perylene molecule resembles 
that of PTCDA. The GGA HOMO-LUMO gap of 1.8 eV 



is slightly larger than that of PTCDA. Fig. 10 gives the 
PDOS on a perylene molecule adsorbed on metal surfaces 
at distances of d — 3.0 and d — 3.6 A. Comparison to 
Fig. |7jb) shows that pinning of the work function sets in 
as the Fermi level reaches the LUMO. In other words, as 
for PTCDA, the work function is pinned by the LUMO 
of the molecule. As the perylene LUMO lies ^ 1 eV 
higher than that of PTCDA, compare Figs. |9] and [lO 




FIG. 10: As Fig. |9] for perylene adsorbed on metal surfaces. 



the pinning level for adsorbed perylene is ~ 3.7 eV. For 
benzene we observe no pinning, see Fig. |6] as the LUMO 
level is too high in energy and is not populated even if 
benzene is adsorbed on a low work function metal like 
Ca. 



B. Charge transfer and interface dipole 

The charge transfer at the PTCDA-mctal interface can 
be visualized directly by calculating the laterally aver- 
aged electron density difference 

An{z) = npTCDA/mctal(2) -"^motal(2:) -^^PTCDA(^)• (9) 

The electron density nmetai and npTCDA of the metal 
substrate and the molecule are obtained in separate cal- 
culations with the substrate and the molecule frozen in 
the adsorption geometry, using the same unit cell. The 
lateral averaging is done as in Eq. (|3|. Examples of An{z) 
are shown in Fig. |11[ 

Fig. ll'a-c) clearly show the formation of interface 
dipoles that are localized at the PTCDA/metal inter- 
face. Note that the sign of the interface dipole moment of 
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FIG. 11: (Color online) Laterally averaged electron density difference An{z) for PTCDA adsorbed on (a) Ca(lll), (b) Al(lll), 
and (c) Au(lll) at a fixed distance d — 3.0 A. The +/— indicate the direction of the interface dipole. The charge on the 
molecule is estimated by integrating over the shaded areas (see text), (d-f) Isodensity surface of An{x,y, z) close to the 
molecular plane. 




PTCDA on Au is opposite to that of PTCDA on the other 
metal surfaces. For PTCDA on Au electrons are dis- 
placed from the molecular region into the metal, whereas 
for PTCDA on other metals electrons are displaced from 
the metal to the molecule. According to the PDOSs 
shown in Figs. |9jc)-(j) the latter can be interpreted as 
electron transfer from the metal to the LUMO of the 
molecule. In contrast, the displacement of electrons for 
PTCDA on Au cannot straightforwardly be related to 
the transfer of electrons to or from a molecular level, as 
can be observed from the PDOSs of Figs. [9]^a),(b). 

The charge displacement for PTCDA on A u can be 
interpreted in terms of the pillow effect .l^^*^ If the in- 
teraction between the molecule and the surface is not 
too strong, then the wave function of the system can be 
written in good approximation as an anti-symmetrized 
product of the wave functions of the separate molecule 
and the substrate. Anti-symmetrization introduces the 
effects of exchange between the molecular and substrate 
electrons, which leads to a repulsive interaction in case 
of a closed shell molecule, the so-called Pauli repulsion. 
The electronic cloud of the molecule is usually "harder" 
than that of the metallic substrate, i.e. the former is less 
easily deformed. The net effect of Pauli repulsion is then 
that the electronic cloud that originally spilled out from 



the metal surface, is pushed back into the metal by the 
molecular electronic cloud (as if the molecule lands on a 
pillow). The push back of electrons into the metal can 
be observed in Fig. [lljc). It lowers the surface dipole, 
and therefore it lowers the work function. 

One would expect that the spatial distribution of the 
displaced electrons reflect the shape of the molecule. 
The pattern of electron depletion in the molecular region 
can be visualized by plotting An{x,y,z), as is shown in 
Fig.[TT|f). The nodal pattern roughly corresponds to the 
second highest occupied molecular vr orbital of PTCDA 
(HOMO-5 in Ref. I117p . It confirms the conclusion drawn 
from the PDOSs of Figs.[9]Ja),(b) that the electron deple- 
tion is not due to a transfer of electrons from the highest 
occupied state of the molecule to the Au surface. 

The change in electron density for PTCDA on other 
metals is qualitatively different from that of PTCDA on 
Au. Electrons are transferred from the metal to states of 
the molecule, which is clearly demonstrated by plotting 
An{x, y, z). The nodal pattern of An{x, y, z) for PTCDA 
on Al, Fig. [iTJe), corresponds to the LUMO of PTCDA, 
indicating that electrons are transferred to this state, in 
agreement with Fig.|9];d). For PTCDA on Ca, Fig.[lT];e), 
the nodal pattern sho ws both features of the LUMO and 
of the LUMO-I- 1.(1221 This agrees with Fig. ^i), which 



13 



TABLE V: Interface dipole per molecule A/i (Eq. (lOl) and 
total molecular charge q (Eq. ( 11 l)for PTCDA on metal sur- 
faces. The distance between the molecules and the surface is 
fixed at 3.0 A. 



metal 



Am (d) 



Au 
Ag 
Al 

Mg 

Ca 



-1.58 
0.77 
2.00 
2.29 
3.29 



+0.34 
-0.31 
-0.82 
-0.84 
-1.34 



general mechanism that should be operative for any ad- 
sorbed closed shell molecule, even if electrons are trans- 
ferred from the substrate to the LUMO, as for PTCDA 
adsorbed on other metal surfaces. One would expect that 
this leads to an interface dipole that depends strongly 
on the molecule-surface distance. However, Fig. [5] shows 
that for PTCDA on Ag and Al (and to a lesser extend 
also for PTCDA on Mg), the work function, and there- 
fore the interface dipole, is independent of the distance 
over a considerable range. 



shows that electrons are transferred to both these states. 

One can calculate the interface dipole per adsorbed 
molecule as A/x = eA J zAn{z)dz, with A the surface 
area of the adsorbed molecule. Alternatively the interface 
dipole is extracted from the change in the work function 
AW upon adsorption of the moleculed^^Sl 



A/i : 



eoA 



AW. 



(10) 



The results are given in Table |V] The total charg e g on 
the molecule can be estimated by 



(11) 



e / An{z)dz, 



where zq is the point between the molecule and the sur- 
face where An(zo) = 0. The integration is indicated by 



shaded areas in Figs, ll'a-c) 



Comparing the values for Ag to Ca in Table [Vj one 
notices an increase in the interface dipole and in the 
number of electrons transferred from the metal surface 
to the molecule. This is consistent with the change in 
the work function upon adsorption and with the PDOSs, 
see Figs. [S] and |8] The charge distributions of the in- 
terface dipole on Ag, Al, Mg, and Ca are similar. The 
distributions shown in Figs. ll^a,b) can be interpreted 
as electronic charge placed in a 7r-orbital on the PTCDA 
molecule, which is screened by electrons in the metal, 
which leads to characteristic Friedel oscillations in the 
metal substrate. 

Adsorption of PTCDA on Au is qualitatively differ- 
ent from adsorption on the other metal surfaces, as can 
be judged from the signs of the interface dipoles and 
the charges. The charge distribution of PTCDA on Au, 



Fig. 11 'c), is also qualitatively different. It is more local- 



ized in the region between the molecule and the surface, 
as can be expected if it is due to Pauli repulsion (i.e. the 
pillow effect), since the latter is active in the region where 
the molecular and the surface wave functions overlap. 

Since this overlap decreases with increasing distance 
between the molecule and the surface, one expects the 
charge displacement to decrease accordingly. Indeed the 
charge on PTCDA adsorbed on Au, calculated using 
Eq. ( 11 1, at a distance d = 3.6 A is 0.23e, as compared to 
0.34e at d = 3.0 A, see Table|vl The pillow effect IS a very 



C. Model 

In this section we aim at setting up a simple model 
that explains qualitatively the most prominent features 
of the work functions of molecular monolayers adsorbed 
on metal surfaces. In particular, we want to identify 
the different regimes, i.e. work function decrease by the 
pillow effect and work function pinning by charge trans- 
fer. Moreover, the lack of the work function dependence 
on the molecule surface distance in the pinning regime 
should be clarified. 

To construct this model we assume that the relevant 
energy scale is set by molecular properties such as the 
HOMO-LUMO gap, and neglect the broadening of molec- 
ular levels introduced by the interaction between the 
molecule and the surface. We assume that electrons can 
be transferred to the LUMO, whose energy eM{N) de- 
pends on the occupation number N , as in Eq. ^. In 
the spirit of Slater's transition state approach we assume 
that N can take any value < < 2.^^ The charging 
energy U depends upon the environment of the molecule. 
All molecules in a monolayer have the same occupation 
number and in the effective U the interactions between 
all molecules should be taken into account, as well as the 
interactions with the metal substrate. 

If the molecular layer is in equilibrium with the metal 
substrate, one has cm^N) = Ep, with Ep the Fermi level 
of the metal. This determines the occupation number 



N = 



Ep — £0 
U ' 



(12) 



The idea is illustrated by Fig 12 ». If iV ^ 0, the molec 



ular layer is charged and with the screening charge in the 
metal substrate this leads to a dipole layer. It results in 
a potential step at the interface, which we parameterize 
as NF with F the potential step normalized per elec- 
tron transferred to a molecule. The work function then 
becomes 



W=Wc + NF, 



(13) 



with Wc = —Ep the work function of the clean metal 
surface. 

So far we have not yet taken the pillow effect into ac- 
count, which lowers the work function of the clean metal 
surface from Wc to Wc — A. As discussed in the previous 
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of Eqs. (14 1 and (15) can be simplified considerably if 



molecule 




LUMO 



HOMO 



molecule 



LUMO 



HOMO 



FIG. 12: (Color online) (a) The pillow effect results in a po- 
tential step A, which lowers the molecular levels close to the 
interface, (b) Charge transfer to the LUMO rises the molec- 
ular levels and pins the Fermi level. 



section, the pillow effect leads to an interface dipole that 
is localized mainly between the surface and the molecule, 
see Fig. 



11 'c). Therefore, it is consistent to apply the 



potential step to the molecular levels and replace eo by 
eo - A in Eq. The idea is illustrated by Fig [ifta) . 

Replacing Wc with W^^.- A in Eq. ([T3| and using Eg. (l2p 
then gives 



W=iWc~A){l 



F 



U' U 

The S parameter, Eq. ([T|), is then given by 



1-^. 
U 



(14) 



(15) 



Since < N < 2, these expressions are valid if A — 
eo — 2U < Wc < A — cq. If the clean metal work function 
falls outside these bounds, one has 



W ^Wr-A: 5* = 1, 



(16) 



for Wc > A — eo- The molecular level is then unoccupied, 
i.e. = and only the pillow effect is operative. If 

Wc < A — eo — 2U, one has N ^ 2. The molecular level 
is fully occupied, which leads to 



W = Wc - A + 2F; S = l. 



(17) 



From purely electrostatic considerations one has < 
F < [/,136 and < S' < 1, cf. Eq. The expressions 



charge on the molecular layer and its counter charge in 
the metal substrate can be modeled by a plane capacitor. 
In that case we have U = F = e^/ C, where C is the 
"capacitance" of the molecule jl^^SEMl -which leads to the 
simple expressions 



W = -eo; 5 = 0. 



(18) 



In this limit one has perfect pinning, i.e. the work func- 
tion is determined by the molecular level only. It is in- 
dependent of the metal and of the pillow effect. 

This simple model can be used to qualitatively describe 
the work functions in Figs. |4]- [7] The simplest case is if 
Wc > A — eo, i.e. if the work functions of all the metals 
considered are too high with respect to the position of 
the LUMO level. W and S are then given by Eq. (16 1, 
i.e. the work function is simply shifted with respect to 
the work function of the clean metal surface. Benzene 
adsorbed on metal surfaces is such a case, see Fig. |6] 
Since the work function shift is determined by the pillow 
effect, one expects it to be sensitive the distance between 
the molecule and the surface, which can be observed in 
Fig.i 

If IVc < A - eo, the LUMO reaches the Fermi level of 
the metal. W and S are given by Eq. ( 14 1 , and in the 



simple plane capacitor model by Eq. ( 18 ) close-packed 



monolayer of planar molecules, such as in the herringbone 
structure of PTCDA (Fig.[ljb)), comes closest to a plane 
capacitor. Fig. [5] shows that indeed the work function is 
pinned at ~ 4.7 eV over a considerable range of metal 
work functions and molecule-surface distances. If Wc < 
A — eo — 2C/, the LUMO becomes fully occupied. W and 
S are given by Eq. (17 1 and the work function becomes 
"unpinned" . 

The energy at which this occurs depends upon the 
molecule-surface distance. Decreasing the distance in- 
creases the pillow effect, i.e. it increases A. Moreover 
it decreases U, since at a shorter distance the screening 
by the metal substrate is larger. In the plane capacitor 
model U oc l/d, where d is the molecule-surface distance. 
The distance dependence of the unpinning of the work 
function is observed in Fig. [H] At d = 3.6 A the work 
function is fully pinned, at d = 3.3 A it becomes un- 
pinned for Ca, and at d = 3.0 A it is unpinned for Ca 
and Mg. 

The dilute structure for PTCDA has only 
coverage (Fig. [ija)). This situation cannot be describe 
by a simple plane capacitor, and one has to use Eq. ( [l4| . 
It gives a linear dependence of W on Wc with a slope 
< S < 1, which can be observed in Fig. |4] 

The behavior of perylene is consistent with the model 
given above. If Wc > A — eg, only the pillow effect is 
operative and the LUMO is unoccupied, cf. Eq. (16 1. 
This holds for Au and Ag in Figs, ^a.) and (b). If Wc < 
A — eo, we observe pinning. For a close-packed monolayer 
the plane capacitor model explains the pinning (Eq. ( |18[ )) 
observed for Mg and Ca in Fig. ^h). For the dilute 
structure with ^ ML packing density, Eq. (14) can be 



i ML 
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used to describe the behavior for adsorption on the low 
work function metals. 



VI. SUMMARY AND CONCLUSIONS 

We study the interface dipole formation at interfaces 
formed by a monolayers of PTCDA, benzene and pery- 
lene molecules with the Au, Ag, Al, and Ca(lll) and 
the Mg(OOOl) surfaces, using first-principles DFT calcu- 
lations. The interface dipoles are monitored by calculat- 
ing the change of the surface work function upon adsorp- 
tion of the molecular layer. Molecular packing densities 
corresponding to | ML and 1 ML coverage are considered 
and the distance between the molecules and the surfaces 
is varied to establish the dependence of the work function 
on these parameters. 

Adsorption of PTCDA in a densely pack structure 
leads to pinning (S = 0) of the work function at ~ 4.7 eV 
for a range of metal substrates and molecule-substrate 
distances, in good agreement with experimental obser- 
vations. The interface dipoles that are created upon ad- 
sorption compensate for the differences between the work 
functions of the different clean metal surfaces. Along the 
series Ag, Al, Mg, Ca the interface dipole generated by 
PTCDA adsorption increases, which is consistent with an 
increasing transfer of electrons from the substrate to the 
PTCDA molecules. The increased electron transfer also 
leads to a stronger bond between the molecule and the 
surface. Decreasing the packing density of the PTCDA 
molecules to | ML decreases the pinning effect, but it 
still gives a linear dependence between the work function 
of the adsorbed layer and that of the clean metal surfaces 
{S « 0.5). Adsorption of PTCDA on Au(lll) leads to 
a very weak bond and an interface dipole that has an 
opposite sign, as compared to the other surfaces. Here 
the pillow effect is dominant, which pushes the electrons 
into the metal substrate. 

Adsorption of benzene results in a reduction of the 



work function, irrespective of the substrate, in agreement 
with experiments. This reduction is in the range 0.2-0.8 
eV, depending upon the distance between the molecule 
and the surface. At a fixed distance S « 0.9. In the case 
of benzene adsorption only the pillow effect is operative. 
The latter, as well as the effect of charge transfer to the 
molecule are observed in adsorption of perylene. Ad- 
sorption of a full ML of perylene molecules on low work 
function metals gives work function pinning (S = 0) at 
~ 3.7 eV. Adsorption on high work function metals gives 
the work function reduction characteristic of the pillow 
effect with S « 0.9. The transition between the two 
regimes takes place for substrate work functions in the 
range Mg-Al. Decreasing the packing density decreases 
the pinning in the low work function regime. 

A simple model inspired by Slater's transition state 
approach allows us to describe the changes in the work 
function upon adsorption qualitatively. The model incor- 
porates the charge transfer between the substrate and the 
molecular levels, the charging energy of the molecules, 
the pillow efi^ect, and the interface dipole layer. It shows 
that for planar molecules the work function is pinned at a 
level that is determined by the molecules and not by the 
substrate. That level does not correspond to the molecu- 
lar EA, but lies within the transport gap of the molecular 
material. 
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